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ABSTRACT 

Identifying the electromagnetic counterparts of gravitational wave (GW) 
sources detected by upcoming networks of advanced ground-based interferom- 
eters will be challenging due in part to the large number of unrelated astrophys- 
ical transients within the ~ 10 — 100 deg 2 sky localizations. A potential way to 
greatly reduce the number of such false positives is to limit detailed follow-up 
to only those candidates near galaxies within the GW sensitivity range of ~ 200 
Mpc for binary neutron star mergers. Such a strategy is currently hindered by 
the fact that galaxy catalogs are grossly incomplete within this volume. Here we 
compare two methods for completing the local galaxy catalog: (1) a narrow-band 
Ha imaging survey; and (2) an H I emission line radio survey. Using Ha fluxes, 
stellar masses (M*), and star formation rates (SFR) from galaxies in the Sloan 
Digital Sky Survey (SDSS), combined with H I data from the GALEX Arecibo 
SDSS Survey and the Herschel Reference Survey, we estimate that a Ha survey 
with a luminosity sensitivity of Lh q = 10 40 erg s -1 at 200 Mpc could achieve a 
completeness of f^ R ~ 75% with respect to total SFR, but only fyg ~ 33% with 
respect to M± (due to lack of sensitivity to early- type galaxies). These numbers 
are significantly lower than those achieved by an idealized spectroscopic survey 
due to the loss of Ha flux resulting from resolving out nearby galaxies and the 
inability to correct for the underlying stellar continuum. An H I survey with 
sensitivity similar to the proposed WALLABY survey on ASKAP could achieve 
/sfr ~ 80% and f^] ~ 50%, somewhat higher than that of the Ha survey. Fi- 
nally, both Ha and H I surveys should achieve > 50% completeness with respect 
to the host galaxies of short duration gamma-ray bursts, which may trace the 
population of binary neutron star mergers. 
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Subject headings: gamma rays: bursts— gravitational waves— galaxies: distances 
and redshifts— radio lines: galaxies 



Introduction 



The inspiral and coalescence of neutron star binaries (NS-NS) are the most likely astro- 
physical sources for direct detection with the upcoming advanced networks of ground-based 



Abramovici et al. 


1992; 


Caron et al. 


1999; 


Acernese et al. 


2009; 


Abadie et al. 



2010) . Maximizing the science achievable from such detections will require t he identifica 



tion and localization of an associated electro magnetic (EM) counterpart (e.g.. iBloom et al. 



20091 : |Phinnevil2009l : Metzger fc Bergeil 12012 . hereafter MB 12). The EM-GW link is crucial 
for identifying the host galaxy and distance of the merger, for placing the mergers in an 
astrophysical context, for studying the hydrodynamics of matter during the merger process, 
and potentially for li f ting degeneracies a ssociated with the inferred binary parameters (e.g., 
Hughes fc Ho3l2003l iKelley et al.ll2012h . 



One commonly-d iscussed EM count erpar t of NS-NS me rgers is a short- duration gamma- 
ray burst (GRB); see iFong et al.l (120101 ) and iBergerl (120111 ) for evidence favoring this asso- 
ciation. Unfortunately, on-axis short GRBs occur at a lo w rate of < 1 yr 1 within the 



200 Mpc range for NS-NS detections by aLIGO/ Virgo (INakar et all 12009: IFong et al. 



2010c MB 12), so the joint detection of a short GRB with a GW event is expected to be rare - 



More isotropic optical counterp arts, such as an off- axis GRB afterglow ( jCoward et al.ll2011 
van Eerten fc MacFadyenl 12011c MB12) or supernova-like emission powered by the decay of 



radioactive e.iecta ( "kilonova" ; iLi fc Paczvhskj|1998c iMetzger et al.ll2010c iRoberts et al.ll2011 



Gorielv et aDboill : iRosswog et ap2012[ ) , may instead represent more promising counterparts 
for the bulk of GW-detected events (MB12); the typical timescale for these counterparts is a 
few days. Delayed radio emission, powered by the interaction of (non-)relativistic ejecta with 
the surrounding environment, may also be detectable if the merger occurs in a sufficiently 



dense medium, with a timescale of months to decades (INakar fc Piranl 12011 
20121 ). 



Piran et al. 



Regardless of physical origin, a major challenge to identifying EM counterp arts is the 
expected poor sky localiza t ions o f ~ 10 — 100 deg 2 for networks of GW detectors (IFairhurst 



2009c iNissanke et al.l 1201 ll . |2012| ). In the optical band, this large area will necessitate deep 
and rapid follow-up with wide-field survey instruments. Still, an even greater challenge may 
be the large number of false positives in these wide fi elds (e.g., background supe rnovae, 
M dwarf flares, shock break-out events, AGN flares; see iKulkarni Kasliwall 12009c MB12; 
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Nissanke et al.ll2012l ). Given the rapid evolution of the predicted optical signal, such false 
positives need to be quickly eliminated so that candidate counterparts could be followed 
up with deeper photometry or spectroscopy before the transient fades. Although rapid 
follow-up is not as essential at radio frequencies, false positives could also be of concern 
in this case (e.g., radio supernovae, tidal disruption events, AGN flares), given the poorer 
angular resolution of wide-field radio survey instruments and our lack of knowledge about 
the t ransient radio sk y on the timescales and depths of the expected merger counterparts 



ve-g- 



Frail et all 120121 ) . 



The contamination from false positives could be reduced by restricting the search vol- 
ume to_Jocatioj]5 neai^ galaxies within the sensitivity range of aLIGO/Virgo (z < 0.046; 



e.g 



Abadie et al.ll2010f ). This approach has the potential to reduce the search area from 
tens of square degr ees to < 0.5 deg 2 , and hen ce to reduce the number of false positives by 



a factor of ~ 10 3 ( iKulkarni fc Kasliwall 120091 ). However, to implement this strategy with 
success, the census of local galaxies must be reasonably accurate and complete. Unfortu- 
nately, thi s is currently not the case: the galaxy catalogs used for recent LIGO/Virgo GW 
follow-up (lAbadie et al.ll2012l ; Evans et al.l 120121 ) are estimated to be only 60% complet e 



(with respect to £>-band luminosity) at 100 Mpc ( iKopparapu et al.ll2008l ; IWhite et al.ll201ll ). 
in which case the true counterpart (most likely at a larger distance) could easily be missed. 
Furthermore, although 5-band completeness is a reason able proxy for the NS -NS merger 



population if the latter traces current star formation (e.g.. iBelczynski et al.ll2002l ). complete- 



ness with respect to stellar mass is also relevant if a sizable fraction of mergers occur in 



early-type galaxies with older stellar populations (e.g.. lO'Shaughnessy et al.ll2008l ). 



In this paper we compare two approaches for completing the local galaxy catalog. The 
first strategy is a nar row-band Ha survey (§2]), an idea that has been discussed previously 



(e.g. 



Rau et al.l 120091 ) but is fleshed out in detail here. The other strategy, explored here for 
the first time, is a wide- field H I survey ([J2])? as is already being planned as a main science 
driver for future wide-field radio arrays. Our main conclusion is that neither Ha nor H I 
surveys of planned sensitivity are sufficient to fully complete the galaxy catalog, especially 
with respect to stellar mass. However, even with respect to star formation completeness, 
an Ha survey suffers from a significant loss of flux (and hence sensitivity and completeness) 
from spatially resolving the disks of the galaxies and from the effects of Balmer absorption 
in the underlying stellar continuum. In £|2]we discuss the merits of combining Ha and H I 



1 Although it is possible that some mergers will oc cur far outside of the host galaxy of their stellar 



progenitors (e.g., Naravan et al. 1992 ^ Kelley et al. 201of). even offset s as large as tens of kpc (as characterize 



Bergerll20ld iFong et alJboid ) would still reduce the required search 



the observed offsets of short GRBs; 
area substantially as compared to a search conducted with no information on sky position. 
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surveys, and explore the relevance of our results in the context of the host galaxies of short 
GRBs Q UID . 



2. Narrow-Band Ha Imaging Survey 

Most studies in the past have assumed that the rate of NS-NS me rgers in the local 
universe traces i?-band luminosity (jPhinneyl Il99ll ; iKopparapu et al.l 120081 ) , as would be ex- 
pected if the rate of mergers is proportional to the current star formation rate (SFR). This 
naturally led to consideration of narrow-band Ha imaging surveys to complete the local 
galaxy catalog since Ha emission is close l y tied to on-going star formation activity (e.g., 



Kennicuttlll992l iGallagher fc Gibson! 1 19941 iGallego et al.lll995l lKennicuttlll998l : iFuiita et al. 



20031 ) and ima ging surveys can efficiently identify an d help characterize galaxies in a given 
redshift range (IDale et al.ll2008l I201CH ; iLy et al.l 1201 f I ). In this section we address the com- 
pleteness of such a survey as a function of depth in terms of star formation and stellar mass, 
and discuss the associated challenges. 

We first note that while SFR is certainly an important quantity, it is also possible that a 
sizable fraction of mergers instead trace total st ellar mass, as determined by the distribution 
of merger times and the sta r formation history ( Belczvnski et al. 2002 , 20061 ; see Figure 1 1 of 
O'Shaughnessy et al.l 12008 ). iLeibler &: Berger find evidence from the host galaxies of 

short GRBs that merge rs may track a combination of SFR and stellar mass, in rough analogy 
with Type la SNe (e.g.. lScannapieco &: Bildstenll2005l ). Given this possibility, in what follows 
we explore completeness with respect to both stellar mass and SFR independently. 

In addressing the issue of Ha completeness we pursue an empirical approach, rather than 
directly making use of derived relations between Ha luminosity and various galaxy properties. 
Our primary source of data are galaxy catalogs with measured Ha fluxes (-Ffia) provided by 
the MP A-JHU emission lin e analysi $ from the Slo an Digital Sky Survey (SDSS) Data Release 



7 (e.g.. IStrauss et al.ll2002l ; lAbazajian et a l. 2009) . Each galaxy is characterized by a s tellar 



mass (determined by photometric m odeling; iKauffmann et al.ll2003l ; ISalim et al.l 120071 ) and 
SFR (determined from emission lines; iBrinchmann et al.ll2004l ). We correct the measured Ha 
fluxes to account for the finite angular size of the spectroscopic fiber relative to the galaxy by 
scaling the fluxes by the difference of the total r-band magnitude (key word r_petro) t o that 
cont ained in the spectroscopic fiber (keyword r_f iber), as described in lHopkins et al.l (120031 ) 



and IBrinchmann et al.l (120041 ) . This assumes uniform Ha surface brightness (see below for 



further discussion) and allows us to study the completeness achieved by an idealized survey 



^http : //www . mpa-garching . mpg . de/SDSS/DR7/ 
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that picks up the entire Ha flux from the galaxy. We also increase the quoted flux errors by 
a factor of 2.47 (as suggested by the MPA-JHU pipeline) and set a conservative threshold 
on the signal-to-noise ratio of S /N > 5 for an Ha detection; values with S /N < 5 are treated 
as upper limits. Finally, we remove all galaxies with "flagged" SFRs (keyword f lagSFR = 
0) or redshifts (keyword z_warning = 1). 

Figure [T] shows our results for the fraction ( "completeness" ) of the total stellar mass 
(red) and SFR (blue) within 200 Mpc as a function of the survey flux deptrfl (-^lim,Ho) for 3J1 
idealized spectroscopic survey. Completeness is calculated as the fraction of the total SFR 
or the total stellar mass in all galaxies with F Hq , > .Fnm,Hcn normalized to the total SFR or 
mass of all galaxies in the redshift range z < 0.046. As the flux limit decreases below a 
few x 10~ 15 erg s _1 cm~ 2 the completeness becomes increasingly uncertain (shown by the 
hatched uncertainty bands), as increasing fractions of galaxies have only upper limits on 
-Fhq. The miminum (maximum) completeness in this case is calculated by assuming that 
none (all) would be detections at Fn a < Fi im; Ha- 

An idealized survey to a depth Fi imi H a = 2 x 10~ 15 erg s -1 cm -2 (a luminosity of 
-^lim.Ha = 10 40 erg s -1 at 200 Mpc), as could be reached with a few minu te integration pe r 
pointing on a meter-class telescope such as the Palomar Transient Factory (ILaw et al.ll2009l ). 
could achieve a completeness of f^ R ~ 97% with respect to SFR, but only f^" « 60 — 80% 
with respect to total stellar mass. It is not surprising that > .f^° given th e known 



correlation between Ha luminosity and SFR (jKennicuttlll998l ; iKewley et al.l 120021 ) . Lower 
completeness with respect to stellar mass is consistent with the fact that « 60% of mass 
in the local universe is in early-type galaxies (jBell et al.l 120031), while they ac count for only 
~ 10% of the Ha luminosity density at z < 0.046 (e.g. iNakamura et al.ll2003l ). 



For comparison in Figure [T] we also plot the completeness with respec t to -B-Band lu 
minosity using data f rom the 11HUGS survey of galaxies at < 11 Mpc (iLee et al.l 12007 
Kennicutt et al.l 120081 ) . and with respect to tot al Ha luminos i ty (using the z = Ha lu- 
minosity function from iGallego et al.l Il995l ; cf. iLy et al.l l2007t iDale et al.l |2010| ). I?-Band 
completeness (f^ a ) tracks closer to SFR than to stellar mass, with fg a f» 87% completeness 
at Ln a = 10 40 erg s _1 . This is not unexpected, as blue light tracks ongoing star formation, 
although not as tightly as Ha luminosity. 

The estimates of mass and SFR completeness in Figure [1] are only appropriate for an 
idealized survey, such as that realized through direct spectroscopy of the target galaxies 



3 We do not discuss here whether the Ha survey should be f ramed as being defined by a flux or equivalent 
width threshold (e.g., see the discussion in Steidel et al. 2000l) . This distinction will have a minor effect on 
our conclusions. 
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(as in SDSS), but with infinite aperture. For an actual imaging surveys, absorption in the 
unde rlying; stellar contin uum can reduce the perceived Ha flux when integrated over a filter 
(e.g., Meurer et all2006). While for young star-forming galaxies this will generally be a minor 



correction (jBrinchmann et al.ll2004[ iMeurer et al.ll2006l ). the correction can be significant for 



early-type galaxies. In Figure [2] we show a comparison between the SDSS spectra of three 
galaxies with similar Ha fluxes, but where the integrated flux in a narrow-band Ha filter 
is reduced by increasing amounts of Balmer absorption from the stellar continuum. In the 
bottom spectrum, the Ha equivalent width (EW) is much greater than the correction due 
to stellar absorption (EW*), in which case the galaxy would be detected regardless of the 
effects of stellar absorption. However, in the top spectrum EW and EW* are comparable, 
and the galaxy might escape detection in a narrow-band imaging survey. 

To explore what completeness could be achieved by an actual imaging survey that 
cannot correct the Ha fluxes for stellar absorption, we recalculated the minimum mass and 
SFR completeness using the fluxes from SDSS which are instead calculated from the line 
equivalent widths computed via integration over broad wavelength ranges. The results are 
shown with dashed lines in Figure [TJ Although the effect on SFR completeness /|fr is 
relatively minor (le ss than a few percent, and comparable to the corrections considered in 



Meurer et al.l 120061 ). the completeness with respect to stellar mass decreases to ~ 47% 



0. 



This difference arises because Balmer absorption is strong and Ha is weak in galaxies with 
relatively old stellar populations and low levels of star formation activity, i.e., those that 
contribute to stellar mass but not SFR completeness. 

Another challenge of a realistic survey is that nearby galaxies generally have angular 
sizes larger than the typical point-spread-function of a ground-based imaging survey. This 
effect can reduce the effective survey sensitivity compared to that for point sources (i.e., 
the survey is actually defined by a limiting surface brightness instead of a limiting flux). 
The actual correction between limiting surface brightness and limiting flux depends on the 
background noise level, the typical seeing, and the distribution of galaxy sizes. In Figure [3] 
we show the distribution of effective angular radii (TZ eS = 7£ ap ,SDSs x 10 (r - f iber ~ r -P etro)/5 ) 
of all galaxies in our sample, where 7£ aPi sDss is the 1.5" radius spectral aperture of SDSS. 
A significant number of sources have radii of > 5", meaning that they will cover > 10 
seeing disks for typical conditions, although this will be somewhat mitigated by the fact 
that Ha typically comes from localized H II regions which have higher-than-average surface 
brightnesses. 

We include the effects of Balmer absorption and the galaxy size distribution to determine 
more realistic completeness fractions. We use the uncorrected SDSS Ha fluxes (i.e., not 
scaling the Ha fluxes by the total r-band magnitude, as described above). This gives an 
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effective seeing disk of 1.5", comparable to typical ground-based observing conditions, and 
basically assumes that the Ha flux outside the central portion is below the detectability 
threshold. We recompute the stellar mass and SFR completeness using this pessimistic 
("imaging") scenario, and the resulting minimum completeness are shown with solid lines 
in Figure [TJ For a (point source) luminosity L^ a = 10 40 erg s -1 , the SFR completeness is 
reduced to f^ R ~ 76%, while the stellar mass completeness is fyg ps 33%. 

In Table[T] we summarize the relevant completeness values. We give three values for the 
Ha; survey, as discussed above. The first corresponds to our ideal survey, where there is no 
Balmer absorption and all of the flux from extended sources is recovered. The other two 
(realistic) sets include Balmer absorption, with the optimistic version still assuming that all 
of the flux from extended sources is recovered, and the pessimistic version assuming that 
none of the flux from extended sources is recovered; the true value likely is between these 
last two versions, i.e. / S H F Q R « 75 - 95% and f$* « 30 - 45% 



3. HI Survey 



Another strategy to complete the local galaxy catalog is via a wide-field H I (21 cm) 
emission line survey, looking for neutral rather than ionized gas. We are motivated in 
particular b y the planned Wide field ASKAP L-band Legacy All-Sky Blind surveY (WAL- 
LABY; e.g., Duffy et al. 2012a ). which plans to observe ~ 75% of the sky0 (declination of 
—90° to +30°) over a timescale of ~ 1 year and detect about a half-million galaxies to 
z ~ 0.26 with an estimated rms sensitivity (over a 100 kHz bandwidth) of ~ 0.7 mJy Be- 
low we examine to what extent WALLABY will be effective at completing the local galaxy 
catalog. Throughout we assume th at H I luminosity (£hi) scales with gas mass (Mm) 
as L H i = 2 x lO^f MHT/Mf^ergs" 1 JDuffv et al] l2012ah. neglecting H I self-absorption (a 
good approximation; IZwaan et al.l 119971 ). Determining the threshold for a given survey also 
requires assumptions about the velocity widths and inclinations of the galaxies. For our 
purposes we assume the limiting flux density and linewidth abov e, which are suffic ient to 
detect M H i ~ 1O 9 M at a distance of 200 Mpc at 5a significance ( IDuffv et al.ll2012aJ Jbh . 



As in the case of Ha, we opt for an empirical approach to estimate the achievable 
completeness with respect to stellar mass and SFR. We are limited by the absence of a 
single large sample of galaxies with measured H I masses (or upper limits), stellar masses, 
and SFRs, yet which is unbiased with respect to galaxy population. For this reason, we use 



4 WALLABY may be complemented by similar surveys in the north, such as the proposed Westerbork 
Northern Sky H I Survey; see http: //www. astron.nl/{~}jozsa/wnshs/survey_layout .html 
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several different H I catalogs to explore completeness with respect to stellar mass and SFR 
(see Figure H] for the H I masses and stellar masses) . 

To study completeness among star forming galaxies, we us e H I masses from the Arecibo 
Legacy Fast ALFA Survey (ALFALFA; iGiovanelli et al.l 120051 ) with SDSS counterparts (for 
which stellar mass and SFR are available as described in §2} in the redshift range z < 0.046. 
Although the entire ALFALFA sample is large (about 9000 galaxies) and covers a wide range 
in stellar mass (M* ~ 10 7 ' 5 — 10 1L5 M ) it is biased towards star forming galaxies and contains 
only SDSS counterparts with H I detections. A complem entary data set is the GALEX 
Arecibo SDSS Survey (GASS; ICatinella et all bold . l2012h . which is much smaller (« 200 
galaxies following our cuts) but also contains deep H I upper limits (down to HI masses of a 
few percent of the stellar mass). Unfortunately, even GASS does not represent a fair sample 
of the entire galaxy population since it is restricted to massive galaxies (M* > 10 10 M M ), 
which for example only account for about half of the short GRB hosts (ILeibler fc Berger 
20101 ); see Figure HI For this re ason, we also use H I masses and upper limits from the 
Herschel Reference Survey (HRS; ICortese et al.ll2012l ). a volume limited sample at < 20 Mpc 
(IBoselli et al.ll2010l ) that covers a relatively broad range in stellar masses (M* ~ 10 8 5 — 10 1L5 
M Q ), covering most of the short GRB host galaxies. 



In Figure [5] we show our results for the completeness achievable by an all-sky H I 
survey with respect to total SFR /gp R (blue) and stellar mass f^ 1 (red), as a function of 
survey limiting flux (-Fi im H i), calculated using the GASS sample (also see Table [I]). For 
mass completeness, we also show the HRS sample (orange). Note that as in the case of Ha, 
completeness is uncertain at low fluxes as the result of H I upper limits, as indicated by the 
hatched regions. 



We find that a survey similar to WALLABY could achieve an SFR completeness of 
/sfr ~ 93 — 98%, but only f^l ~ 44 — 50% with respect to total stellar mass. For mass 
completeness we quote numbers from HRS instead of GASS, since HRS samples a wider range 
of galaxy masses (Fig. H]) . Again, a high SFR comple teness is expected given the correlation 
between cold gas mass and SFR (IBigiel et al.l 120111 ). Lower completeness with respect to 



stellar mass is expected sin ce early-type galaxies are typically deficient in H I (ISage et al. 
2007 ; lOosterloo et al.ll2010l ) yet they contain w 60% of mass in the local universe (jBell et al. 
20031 ). For comparison, we also shows completenes s with respect to total H I mass (as 



determined from the local H I luminosity function of IZwaan et al.l 120051 ) . which is found to 
lie between that of stellar mass and SFR (Fig. [5j Tabled]). 

As we discussed in the context of Ha surveys, for H I surveys we should also consider the 
possible detrimental effect of spatially resolving the galaxies. However, the poorer angular 
resolution of ASKAP (10 — 30") makes this less of an issue than for an optical survey, and 
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the ex pected effect on detectability is anticipated to be minimal (see Figure 5 of iDuffy et al. 



2012a|). The limited angular resolution of the H I surveys will not affect the follow-up of 



GW tri KK crs, £is db ±30" position would be quite sufficient and the H I centroid will be 
known to even better precision. Ideally, candidate galaxies detected by WALLABY will 
also be followed up to obtain optical counterparts b y wide-field sur vey instruments such a s 
VISTA Hemisp heric Survey flAmaboldi et ali l2007r). VST ATLAS f Capaccioli et~al]l2005h . 
or SkyMapper feeller et alil2007h . 



4. Discussion 

Table [1] compares the minimum completeness to SFR and stellar mass achieved by our 
fiducial Ha and H I surveys based on Figures [I] and [5j Overall we find that when considering 
realistic surveys including the effects of Balmer absorption and finite source sizes, an H I 
survey can achieve a somewhat better completeness in both SFR and M* compared to Ha. 

It is clear from our results that neither the proposed Ha or H I surveys will produce a 
local galaxy catalog that is entirely complete with respect to stellar mass out to 200 Mpc. 
Fundamentally, both Ha and H I trace ga s, while ~ 60% of stellar mass in the local universe 



is in early- type galaxies (IBell et al.ll2003l ) with little gas. However, even completeness with 
respect to SFR is unlikely to reach 100% from either survey due to the issues discussed in 
Pand P 

A remaining question is whether there are any advantages to combining the results of 
separate Ha and H I surveys. In Figure E] we show Lh« versus Mhi from the ALFALFA- 
SDSS sample of star forming galaxies. We note that an obvious correlation exists between 
Mhi and Lh q since both trace star formation, but there is also a large scatter (« 1 dex) 
in the correlation where the sensitivity threshold of WALLABY (Mhi ~ 10 9 M Q at 200 
Mpc) intersects that of our fiducial Ha survey (L Hq ~ 10 40 erg s' 1 at 200 Mpc). This 
scatter implies that some fraction of galaxies which are just below the detection threshold 
in H I will be detectable in Ha and vice versa. Therefore the SFR completeness attained 
by combining the results of both Ha and H I surveys together should be more than that of 
either individually. 

If one were to combine H I and Ha surveys, which individually achieve a similar com- 
pleteness (as is the case for WALLABY and our 'ideal' spectroscopic Ha survey)@, then 



5 The fact that a smaller fraction of sources lie in the top-left quadrant (detected in Ha but not H I) 
of Figure |B] than in the bottom-right (detected in H I but not Ha) appears at odds with the comparable 
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from the SDSS- ALFALFA sample we estimate that /sfr could be increased from ~ 95% to 
~ 98%. Mass completeness could be increased by a somewhat greater amount (changing by 
up to ~ 10 — 20%), however _f Mt is ultimately limited to a value < 60% due to the fraction 
of early-type galaxies which possess neither detectable H I or Ha emission. We conclude 
that combining surveys provides the greatest benefit if both H I and Ha surveys separately 
achieve a similar completeness, however even in this case the gains will be relatively modest. 



4.1. Connection to Short GRB Hosts 

So far we have focused on completeness with respect to SFR and stellar mass separately. 
However, it is also of interest to investigate how complete the Ha or H I surveys will be with 
respect to the known population of short GRB host galaxies. We attempt to answer this 
question empirically by investigating the fraction of galaxies with properties (SFR or stellar 
mass M*) similar to those of the short GRB hosts that would be detected in Ha or H I. 



In our calculation we use short GRB 



lost 



and SFRs (or upper limits) from iBergerl (120091 ). resulting in a sample of 11 galaxies. We 



galaxy masses from lLeibler fc Bergerl (120101 ) 



define the 'short GRB host completeness' as the fraction of galaxies in our sample that have 
masses (and in some cases SFRs) within a factor of ± 0.5 dex of those of each short GRB 
host and that are detectable at a given Ha or H I flux, the result of which is then averaged 
over the short GRB hosts. In the case of Ha, the galaxy masses, SFRs, and Ha fluxes are 
again taken from the SDSS sample O, the latter of which are not corrected for stellar 
absorption or galaxy size (as would characterize a purely imaging survey). In the case of H I 
we use galaxy masses from the HRS sample, but we cannot make a cut on the SFR since 
these data are not available. 

Figure [7| shows our results for completeness as a function of Ha or H I flux, normalized 
to the sensitivity of our fiducial surveys. In the case of Ha, results are shown for two cases: 
(1) in which the galaxy samples are chosen to match both the masses and SFRs of the short 
GRB hosts, and (2) in a case for which the samples are chosen based just on sharing similar 
stellar masses with the short GRB hosts. By this criteria, we find that our fiducial Ha 
imaging survey could achieve a short GRB host completeness of /|grb ~ 50% and ~ 25% in 
cases (1) and (2), respectively. Case (1) is more realistic and results in a higher completeness 
because most of the short GRB hosts (9 of 11) are star- forming galaxies, so the SFR cut 
preferentially picks out Ha-luminous galaxies. In both cases the minimum completeness 



completeness attained by our fiducial H I and spectroscopic Ha surveys (cf. Fig.HJEJ; however, the ALFALFA 
sample includes only H I detections and hence overestimates the fraction of H I-detected galaxies. 
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asymptotes to « 80% at low fluxes since 2 of 11 of the short GRB hosts are early-type 
galaxies, for which there are only upper limits on their SFRs and Ha luminosities. 

In the case of H I, we find that our fiducial survey similar to WALLABY could achieve a 
completeness of /|qrb ~ 45%. However, this probably underestimates the true completeness 
for the same reason as with Ha: most of short GRB hos ts are star- forming and there is a 



correlation between SFR and H I mass ( iBigiel et al.ll201ll ). yet no cut was made on SFR 



We conclude that both Ha and H I surveys could achieve at least ~ 50% completion with 
respect to GRB host galaxies, although this number is likely to be substantially higher in the 
case of an H I survey. There is room to increase the completeness by performing a deeper 
search than our fiducial survey, potentially approaching the ~ 80% maximum completeness 
expected if one could detect all 9 of 11 star-forming short GRB hosts at ~ 200 Mpc. 



5. Conclusions 

We used observed samples of nearby galaxies, drawn from SDSS and augmented by H I 
data from several surveys, to estimate the completeness of fiducial Ha and H I surveys to SFR 
and stellar mass out to a distance of 200 Mpc. We conclude that neither H I nor Ha surveys 
to proposed depths will be entirely effective at completing the local galaxy catalog with 
respect to stellar mass within 200 Mpc, but that we can expect reasonable completenesses of 

/sfr ~ 76 %> fwu ~ 33% > /sfr ~ 93 %> and /m! ~ 44 %- At this P° int neither the large-scale 
Ha nor H I surveys discussed here are underway, and it may be possible to alter the strategy 
(i.e., survey depth) to increase the completeness. For instance, halving the survey threshold 
of the fiducial Ha imaging survey results in a substantial increase in the completenesses to 
/sfr ~ 90%, fyg ~ 49%, but for H I only results in a modest increase to /gp R « 95%, 
and /g 1 « 50%. For Ha, these increases are comparable to that achieved by pursuing a 
spectroscopic survey (Fig. [p. Such a deeper survey would also increase the completeness 
with respect to short GRB host galaxies by > 10% (Fig. [7]). 

If the survey thresholds remain as discussed here, assessing whether or not this is a 
concern for aLIGO/ Virgo follow-up largely depends on the fraction of mergers that occur in 
early-type galaxies. If the host galaxies of short GRBs indeed represent a faithful sampling 
of the merger p opulation, th en the current census of about a 5 to 1 ratio of star forming to 



elliptical hosts ( lBergerll201ll ) would suggest that incompleteness of early-type galaxies is not 
a major concern. On the other hand, it is also possible that short GRBs do not represent 
all mergers, or that the short GRB host population is biased against early-type galaxies. 
Current population synthesis models do allow for a sizable fraction > 20 — 50% of elliptical 



hosts at z = ( jO'Shaughnessy et al 



20081 ). 



We note that our treatment of the surveys used simply-defined empirical data sets, 
and some secondary effects may change the results slightly. For instance, neither the Ha 
flux nor the Balmer absorption is expected to be distributed uniformly within the galaxies. 
While we do not think that effects such as these will greatly change our conclusions, future 
studies could be done with more carefully constructed samples or could make use of observed 
semi-empirical correlations between galaxy properties to directly calculate the completeness. 



We thank Jarle Brinchmann and Christy Tremonti for helpful conversations and guid- 
ance using the MPA-SDSS galaxy data; Barbara Cantinella and Luca Cortese for helpful 
conversations and guidance using the GASS and HRS H I data; and Dawn Erb, Wen-fai 
Fong, Alicia Soderberg, and Mansi Kasliwal for helpful discussions. We are grateful to the 
KITP in Santa Barbara for hosting the program "Chirps, Mergers and Explosions", where 
this work began. BDM was supported in part by NASA through Einstein Postdoctoral Fel- 
lowship grant number PF-00065. BDM also acknowledges support from the Lyman Spitzer, 
Jr. Fellowship awarded by the Department of Astrophysical Sciences at Princeton Univer- 
sity. This research was also supported in part by the National Science Foundation under 
Grant Nos. PHY11-25915 (KITP), AST-1008353 (DLK), and AST-1107973 (EB). Funding 
for the SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, the Partic- 
ipating Institutions, the National Science Foundation, the U.S. Department of Energy, the 
National Aeronautics and Space Administration, the Japanese Monbukagakusho, the Max 
Planck Society, and the Higher Education Funding Council for England. The SDSS Web Site 
is http://www.sdss.org/ The SDSS is managed by the Astrophysical Research Consortium 



for the Participating Institutions. The Participating Institutions are the American Museum 
of Natural History, Astrophysical Institute Potsdam, University of Basel, University of Cam- 
bridge, Case Western Reserve University, University of Chicago, Drexel University, Fermilab, 
the Institute for Advanced Study, the Japan Participation Group, Johns Hopkins University, 
the Joint Institute for Nuclear Astrophysics, the Kavli Institute for Particle Astrophysics and 
Cosmology, the Korean Scientist Group, the Chinese Academy of Sciences (LAMOST), Los 
Alamos National Laboratory, the Max-Planck-Institute for Astronomy (MPIA), the Max- 
Planck-Institute for Astrophysics (MPA), New Mexico State University, Ohio State Univer- 
sity, University of Pittsburgh, University of Portsmouth, Princeton University, the United 
States Naval Observatory, and the University of Washington. 



13 




Fig. 1. — Estimated completeness of an all-sky, narrow-band Ha survey with respect to 
total stellar mass (red) and total star formation (blue) as a function of survey depth, Fn mt n a , 
calculated using Ha fluxes and galaxy properties derived from SDSS (see text). Dotted lines 
show the completeness of an idealized spectroscopic survey which measures the entire Ha 
luminosity of the galaxy (infinite aperture) and corrects the Ha flux for stellar absorption; 
the cross-hatched region represents the uncertainties due to Ha upper limits. Dashed lines 
show how the minimum completeness decreases when one does not correct Ha fluxes for the 
underlying stellar Balmer continuum, as appropriate for narrow-band imaging. Solid lines 
show the minimum completeness when the Ha fluxes are also not corrected for the finite 
angular size of the galaxy (assuming a 1.5" radius aperture). These last two cases likely 
bracket the completeness provided by a purely imaging survey. Also shown for comparison 
is compl eteness with respect to S-band luminosity (green) of the local (< UMpc) 11HUGS 



survey (|Kennicutt et al. 



20081 ) and with respect to total Ha luminosity (gray; using the 



Gallego et al.l 1 19951 luminosity function). 
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Fig. 2. — Spectra of three galaxies from our SDSS sample which illustrate, from bottom to 
top, the effect of increasing stellar absorption on the measured Ha equivalent width (EW) 
in a narrow-band survey. The line centers of the Ha, [N II], and [S II] lines are shown 
with red vertical dashed lines (wavelengths are shifted to the rest frame). The gray shaded 
region represents the approximate bandwith of a z = — 0.01 Ha filter. For each galaxy 
we give the EW of the Ha emission line and the EW of the Balmer absorption (EW*), with 
positive values indicating emission; all values have been scaled to have the same continuum. 



-1 —2 

cm , 



58 x 



The Ha fluxes (from the MPA-JHU emission line analysis) are 72 x 10 ergs 
10~ 17 erg s" 1 cm™ 2 , and 37 x 10" 17 erg s" 1 cm" 2 , from bottom to top. It is clear that the 
Ha EW in the bottom galaxy greatly exceeds the Balmer absorption (by a factor of > 30), 
while in the upper galaxies the Balmer absorption is 30% and 100% of the Ha emission, 
even as the line flux changes by only a factor of 2. This would make detecting the top 
galaxy in a narrow-band survey difficult, although the [N II] doublet will generally provide 
a minimu m flux of 10% of the original Ha flux (dependent on metallicity, among other 



factors; IPettini fc Pagelll2004f ). 
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Fig. 3. — Normalized distribution df/d72. e g- of effective angular radii 7£ e g- of galaxies in our 
SDSS sample at redshift z < 0.046 (containing 103,025 sources). The effective radius is 
denned as TZ eS = ft ap , S DSS x 10 (r - f iber - r -P etro )/5 ; wne re r_fiber and r_petro are the 
fiber and Petrosian r-band magnitudes and 7£ aPi sDss is the 1.5" spectral aperture of SDSS. 
To the extent that Ha and r-band have the same surface brightness distribution, the flux 
sensitivity to an extended Ha source is reduced by a factor ~ (7£ e fr/7^app) 2 compared to a 
point source. 
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Fig. 4. — Stellar mass M* versus H I mass Mhi from samples of galaxies used in estimating 
H I completeness: ALFALFA (light blue points, detections only); HRS (brown, both detec- 
tions [squares] and upper limits); and GASS (red, both detections [stars] and upper limits); 
see t ext for references . Vert ical blue lines show the masses of short GRB host galaxies 
from iLeibler fc Bergerl (120101 ). while the horizontal purple line shows the approximate H I 
sensitivity threshold of WALLABY at 200 Mpc. 
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Table 1. Galaxy Completeness within 200 Mpc 



Survey Type 


Proposed Depth 




min(/ SF R) 


min(/ M J 


Ha (ideal spectroscopic 11 ) 


2 x 10~ 15 erg cur 2 s" 


-i 


97% 


58% 




1 x 10 -15 erg cm -2 s" 


-l 


97% 


59% 


Ha (realistic imaging, optimistic 13 ) 


2 x 10~ 15 erg cm~ 2 s~ 


-l 


95% 


47% 




1 x 10~ 15 erg cm~ 2 s~ 


-i 


96% 


52% 


Ha (realistic imaging, pessimistic ) 


2 x 10~ 15 erg cm~ 2 s~ 


-i 


76% 


33% 




1 x 10~ 15 erg cm -2 s~ 


-i 


90% 


49% 


H I 


0.7 mJy x 100 kHz 




93% 


44% 




0.35 mJy x 100 kHz 




95% 


49% 



Note. - We give the minimum completeness to SFR and stellar mass based on our 
strawman Ha and H I surveys. The Ha flux limit corresponds to a luminosity of L-^ a = 
10 40 erg s _1 at 200 Mpc. The H I flux limit corresponds to a luminosity of Lhi = 3 x 10 31 
erg s , or an H I mass of Mhi ~ 10 9 M . 

a This scenario assumes no Balmer absorption and that all of the flux from extended 
sources is recovered (dotted line in Fig. [[]). 

b This scenario includes the effects of Balmer absorption, but assumes that all of the 
flux from extended sources is recovered (dashed line in Fig. [1]). 

c This scenario includes the effects of Balmer absorption, but assumes that none of the 
flux from extended sources is recovered (solid line in Fig. [TJ) . 
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Fig. 5. — Estimated completeness of an all-sky H I survey with respect to total stellar mass 
(red) and total star formation (blue) as a function of survey depth (F^hi), calculated using 
H I fluxes from the GASS sample and galaxy properties derived from the SDSS catalog (see 
Fig. H]text). Survey depth is calculated as the 5a rms sensitivity over an assumed 100 kHz 
bandpass (see text). We also show the mass completeness calculated using H I masses from 
the Herschel Reference Survey (HRS; orange) and the completeness wit h respect to t o tal H I 
mass (brown) calculated using the local H I luminosity function from IZwaan et al.l ( 120051 ). 
The hatched regions indicate the range of uncertainty in completeness at low flux due to H I 
upper limits. 
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Fig. 6. — Ha luminosities L^ a versus H I masses Mm from the ALFALFA-SDSS sample 
of star- forming galaxies, showing both Ha detections {blue stars) and upper limits {red 
triangles). Although a clear correlation exists between L^ a and Mhi (since both trace star 
formation), the scatter in Lu Q is significant (~ 1 dex) near the sensitivity threshold of an H I 
survey similar to WALLABY {purple line). The SFR completeness obtained by combining 
Ha and H I surveys may thus be increased somewhat compared to that obtained by either 
individually. 
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Fig. 7. — Completeness of Ha and HI surveys with respect to galaxies with properties (SFR 
and M*) similar to those of the host galaxies of short GRBs, normalized to the sensitivity 
of our fiducial Ha imaging survey and to that of WALLABY. In the case of Ha, the SDSS 
galaxy sample is used; results are shown both using a subsample of galaxies selected based 
on similar stellar masses and SFRs to the short GRB hosts (brown), as well as subsamples 
chosen based just on similar stellar masses (blue). In the case of H I, the HRS sample is used 
and results are shown just for the subsample with similar stellar masses (orange); since most 
SGRB hosts are star forming, the completeness achievable by H I is probably underestimated 
by this figure (see text). 
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